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ABSTRACT 


The  problem  of  dynamic  stresses  and  displacements  around  a  cavity 
and  rigid  inclusion  of  arbitrary  density  is  examined  for  an  elastic 
medium  during  the  passage  of  a  plane  shear  wave.  In  the  cavity  case, 
the  dynamic  stresses  and  displacements  are  found  to  be  dependent  upon 
the  incident  wave  number  and  Poisson’s  ratio  of  the  medium.  In  the 
rigid  inclusion  case,  it  is  found  that  dynamic  stresses  and  the  rigid 
body  rotation  and  translation  are  dependent  upon  the  incident  wave 
numbers,  the  Poisson’s  ratio,  and  the  density  ratio  of  the  medium  and 
the  insert.  Close  coupling  is  observed  between  the  stresses  and  the 
rigid  body  motion  of  the  insert. 
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INTRODUCTION 

The  problem  of  dynamic  stress  concentration,  around  a  cylindrical 

discontinuity  in  an  infinite  elastic  medium,  subjected  to  a  plane  travelling 

r  1  2  3  5] 

wave  has  been  under  extensive  study  in  recent  years.  >  >  >  >  J 

Former  investigations,  however,  considered  only  the  case  where  the  incident 
wave  is  a  plane  compressional  wave.  In  this  paper,  the  dynamic  stresses 
due  to  an  incident,  plane,  shear  wave  are  investigated.  Two  cylindrical 
discontinuities  are  considered:  a  cavity,  and  a  rigid  inclusion  of 
arbitrary  density.  Furthermore,  the  displacements  on  the  cavity's  free 
surface  and  the  rigid  body  translation  and  rotation  of  the  rigid  inclusion 
are  also  determined. 

In  the  approach  presented  here,  the  incident  shear  wave  is  assumed 
to  vary  harmonically  with  time.  At  the  boundary  of  the  discontinuity, 
incident  waves  will  be  reflected  into  compressional  and  shear  waves.  The 
resultant  dynamic  stresses  around  the  cylindrical  discontinuity  can  be 
found  frbm  the  solutions  of  the  scattered  waves  and  the  incident  waves 
near  the  obstacle.  The  solutions  obtained  are  based  on  the  generalized 
plane  strain  assumption. 

Comparisons  are  made  between  Kirsch's  static  solution  for  the 
cavity  case  and  Goodier's  ^ ^  static  solution  for  the  rigid  inclusion 
with  the  present  solution.  As  the  wave  numbers  become  small,  the  dynamic 
cavity  solution  reduces  to  the  static  solution.  The  dynamic  solution  for 
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the  rigid  inclusion  also  reduces  to  the  static  solution  for  any  finite 

inclusion  density.  However,  when  the  inclusion  density  is  infinite 

(which  is  the  same  as  holding  the  inclusion  fixed  in  space)  a  singularity 

will  occur  as  the  wave  numbers  approach  zero.  A  similar  phenomenon  was 

[4] 

discussed  in  detail  by  Pao  and  Mow. 


GENERAL  THEORY 

Under  the  generalized  plane  strain  assumption,  the  displacement 
equations  of  motion  are 

(X  +  p)  V  V  .  u  +  [i  u  .  p  u  (1) 


where 

u  is  ,the  displacement  vector 

V  is  the  vector  differential  operator 
X  and  fju  are  the  Lame  constants 
p  is  the  density 


The  displacement  vector  u,  under  the  Helmholtz  theorm,  can  be 
represented  in  terms  of  a  scalar  potential  and  a  vector  potential.  In 
the  case  of  plane  strain,  this  displacement  vector  is  given  by 

u  -  V  cp  +  V  X  (e  i|r) 


(2) 


3 


where  e  is  a  unit  vector  along  the  axis  of  the  cylindrical  discon¬ 
tinuity.  Each  potential  function  then  satisfies  a  scalar  wave  equation 


a 


2 

V  9  =  9 


(3) 


C 


2 

3 


v  ♦  -  ♦ 


(4) 


where 


c2  =  k  ±  2U 
a  p 


INCIDENT  AND  REFLECTED  WAVES 

An  incident  plane  shear  wave  propagating  in  the  positive  x 
direction  is  represented  by 


*(i)  =  t0e1(Px  "  (Ut) 

«p<4>  "  0 


(5) 


k 


where 

^  is  the  measure  of  amplitude 
oj  is  the  circular  frequency 

3  =  <u/c  is  the  wave  number  of  the  shear  wave 

P 

Equation  (5)  can  be  expressed  in  polar  coordinates 

00 

f(i)  =  tn  y  e  inJ  Or)  cos  n0e“iu)t  (6) 

u  n  n 

n=0 


where 


1  n=0 
.2  rt>l 


If  the  origin  of  the  polar  coordinates  is  selected  at  the  center 
of  the  cylindrical  discontinuity  (Fig.  1),  the  reflected  waves  can 
be  expressed  as 


cp(r)  =  J  An  Hn  s*n  n®e 


-iu)t 


n=0 


=  Y  Bn  Hn  cos  nee‘lu)t 

n=0 


(7) 


(8) 
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where 


A  and  B  are  expansion  coefficients  to  be  determined  by  the 
n  n 

appropriate  boundary  conditions 


a  =  (jd/c  is  the  compress ional  wave  number 
a 

H  denotes  the  Hankel  function  of  the  first  kind  of  order  n. 
n 


CYLINDRICAL  CAVITY 


RIGID  CYLINDRICAL  INCLUSION 


Fig.  1 


The  total  potential  can  then  be  determined  by  superimposing  the 
incident  and  reflected  potentials.  Thus,  the  displacement  potentials 
are 


(i)  A  (r)  (r) 

cp  =  cp  +  cp  =  cp 


,  ,(i)  ^  , (r) 


(9) 


Substitution  of  Equations  (6)  -  (9)  into  the  well-known  relation* 

fAl 

ships  between  stresses,  displacements  and  displacement  potentials 
yields  the  corresponding  stress  and  displacement  components .  With  the 
time  factor  e  omitted,  they  are  given  by 
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rr 


-  2\ir  2  y  (i|fne  innD  -I-  A  D  +  B  0D  )  sin  n0 
U  T0  n  0  In  n  1  In  n  2  In 


n=0 


T00  =  2^r’2  1  (^Oenin  0°2n  +  An  lD2n  +  Bn  2°2n)  8in  n® 


n=0 


(10) 


r  =  2p,r  2  Y  (iji  e  in  -.D  +  A  D  +  B  „D  )  cos  n0 

r0  L,  T0  n  0  3n  n  1  3n  n  2  3n 


n=0 


and 


-r  1  V  (fAe  in  nD,  -  A  D.  +  B  0D,  )  sin  n0 
Zj  On  0  4n  n  1  4n  n  2  4n 


n=0 


-r"1  y  (i|rne  in  -.D.  -  A  ,0.  +  B  ,D_  )  cos  n0 

Zj  y0  n  0  5n  n  1  5n  n  2  5n 


n=0 


(ID 


In  Equations  (10)  and  (11),  AD;  through  AD_  represents 

u  in  u  on 

the  contribution  on  stresses  and  displacements  by  the  incident  wave, 
while  ~  2D5n  and  2°ln  "  2°5n  rePresents  the  contribution  on 

stresses  and  displacements  by  the  reflected  compress ional  and  shear 
waves,  respectively.  These  functions  are  defined 
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0Dln  =  n^n+1)  Jn  "  nPr  Jn-i 

0°2n  =  Jn  +  nPr  Jn_i  (3r) 

2  2 

0D3n  =  '(n2+n  "  V"  }  Jn(Pr)  +  er  Jn-1  (3r) 

0D4n  "  nJn  <**> 

0D5n  =  -nJn  Or)  +  Pr  Or) 

2  2 

-D  =  (n2+n  -  )  H  (ar)-ar  H  (ar) 

X  in  in  n-1 

2  2  2  a2  2 

^2n  =  “(n  +n  -a  r  +  (ar)  +  ar  H  ,  (ar) 

jD3n  =  -(n2+n)  Hn  (ar)  +  nar  H  j  (ar) 

lD4n  =  -n  Hn  (ar)  +  ar  Hn-1  (ar) 

,D  =  n  H  (ar) 

1  5n  n 


(12) 


(13) 
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2Dln  =  n(^nfl)  Hn  "  n^r  Hn_i  Or) 

2D2n  =  -n (n+1)  Hn  (gr)  +  ngr  (gr) 

2  2 

2D3n  =  -(n2+n  -  (gr)  +  gr  (gr)  (14) 

2°4n  ‘  ®'> 

2D5n  ‘  -*n  ««  *  Sc  Vl  ®r) 


BOUNDARY  CONDITIONS  AND  SOLUTIONS 

As  previously  mentioned,  the  coefficients  A  and  B  are  to 

n  n 

be  determined  by  the  appropriate  boundary  conditions  associated  with 

the  two  types  of  discontinuities  considered:  the  cavity  and  the  rigid 

inclusion.  The  appropriate  boundary  condition  for  the  cavity  case  is 

obvious;  i.e.,  traction  free  at  the  boundary.  The  boundary  condition 

[41 

for  the  rigid  inclusion,  however,  is  not  this  clear.  The  correct 
boundary  condition  in  this  case  is  obtained  by  leaving  the  rigid  in¬ 
clusion  free  to  rotate  as  well  as  translate  with  the  surrounding  medium. 
This  implies  that,  at  the  interface,  the  displacements  due  to  the  combined 
incident  and  reflected  waves  are  equal  to  the  displacements  corresponding 
to  the  rigid  body  rotation  and  translation  of  the  inclusion  which  are 
determined  by  the  application  of  Newton’s  law  of  motion. 
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Cavity 

The  boundary  conditions  at  r  *  a,  where  a  is  the  radius  of  the 
cylindrical  discontinuity,  are  given  by 


T 

rr 


r@ 


=  0 


If  each  term  of  the  series  for  t  and  T  „  in  Equation  (10)  is 

rr  r0 

set  equal  to  zero  at  r  *  a,  the  coefficients  A  and  B  are 
n  1  n  n 

determined  from  a  pair  of  equations  for  each  value  of  n 


A 

n 


Vn1 


(An  26ln 
(An  263n 


(15) 


B 

n 


Vn1 


An  (An 
An  (An 


and 


A  * 
n 


lDln  2°ln 


l°3n  2D3n 


(16) 


(17) 
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where  J)  . denote  the  value  of  . „D. 

0  In  2  3n  0  In  2  3n 

at  r  =  a.  The  displacements  and  stress  can  then  be  calculated  by  the  sub¬ 
stitution  of  Equations  (15)  through  (17)  into  Equations  (10)  and  (11), 
and  the  solution  is  complete. 


It  is  convenient  to  nondimens ionalize  the  stresses  and  dis¬ 
placements  with  respect  to  the  dimensional  coefficient  of  the  stress 

2 

and  displacement  associated  with  the  incident  wave,  i.e.,  =  p,0  i|tq 

and  uQ  =  -ip^o  *  *s  not®^  that  the  wave  number  is  normalized  by 

the  radius  of  the  circular  inclusion.  Thus,  the  dimensionless  wave 

numbers  aa  and  |3a  may  be  regarded  as  the  ratio  of  the  circumference 

to  the  incident  wave  length.  Thus,  at  the  boundary  (r  *=  a)  the 

*  *  * 

dimensionless  stress  g  and  displacements  u^  and  Ug  are  given  in 
the  following  form 


ee 


-5«-W  ‘ 


n(„2-l.  ^  >  H  fa.) 


sin  n0 


n=l 


(18) 


*  2Ba  r  nin  [ (n+l)Hn(oa)-  oa  Hn_1(aa)] 
u  =  -  — >  - 7 - 

r  n  L  A 

1  A 

n=l 


sin  nd 


(19) 
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e  in  [ (n^+n  -  ®9a  )  H  (aa)-  aa  H  (aa) 

-S— - — a - n=I -  cos  ne 


(20) 


where 


111-21 
(1-2  a) 


a  is  the  Poisson’s  ratio. 


_  jL 

2 

a 


a 


00 


n=0 


Rigid  Inclusion  of  Arbitrary  Density 

As  previously  discussed,  the  rigid  inclusion  in  general  will 
rotate  and  translate  with  the  surrounding  medium.  The  resultant  forces 
and  torque  acting  on  the  inclusion  as  a  result  of  the  stresses  acting 
on  the  boundary  due  to  the  incident  and  reflected  waves  are  expressed 
by 


2h 

F  =  I*  (t  cos  0  -  t  a  sin  0)  a  d0 
x  Jq  rr  r0 


2tt 


F  =  f  (t  sin  0  +  t  A  cos  0)  a  d0 
v  v  ^  rr  r9 


=  / 


2n 


t  &  a  d0 

re 


(21) 


(22) 


T 


(23) 
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By  substituting  Equation  (10)  with  r  «  a  into  Equations  (21)  through 
(23)  and  considering  the  orthogonality,  it  can  be  readily  shown  that 


F  -  0 

x 


(24) 


F 

y 


2[ia 


"  1  <0511  +  0531>  +  *1 

+  B1  ^2B11  +  2B31^ 


<1D11 


*  iV 


(25) 


and 


T  =  4rv  Ct0  0D30  +  Aq  1D30  + 


B0  2D30-^ 


From  Newton's  law  of  motion,  it  is  known  that 


(na2p1)  U 


F 

y 


C^) 


T 


(26) 


(27) 

(28) 


By  restoring  the  time  factor  e  which  is  omitted  in  Equation  (10), 

and  twice  integrating  Equations  (27)  and  (28)  with  respect  to  time, 
the  rigid  body  translation  and  rotation  of  the  inclusion  is  found  to 
be  given  by 
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uy  ■  r  -fr  t  2V  <05U  +  0531>  +  A1  <15U  +  1531> 

p  a 


+  B1  (2511  +  2631)  ^ 


(29) 


0  ■  -Jib  %  oSo  +  A0  1°30  +  B0  2°30) 

P  a 


(30) 


where  T\  = 


is  the  ratio  of  density  of  the  surrounding  medium 


to  that  of  the  inclusion.  The  translation  and  the  rotation  will 
vanish  if  T|  =*  0,  Hence,  the  problem  is  the  same  as  the  case  when  the 
inclusion  is  fixed  in  space. 

For  any  inclusion  with  finite  density,  the  displacements  at  the 
boundary  of  the  inclusion  must  be 


u  =  U  sin  0 

r  y 


(31) 


=  U  cos  0  +  a© 
0  y 


If  the  expressions  for  u^  and  u^  in  Equation  (11)  are  substituted 
into  Equation  (31)  and  the  orthogonality  conditions  are  employed,  it  is 


found  that  for 


Ik 


e2a2  J.(pa)  +  8T|  [  &  J  (Pa)  -  J.  (pa)  ] 

- i - i - H - i -  (32) 

pV  H^pa)  +  8^  [  ^  HQ(Pa)  -  H^pa)  ] 


n  -  1 

2i'tlo 

Ax  - - -  pa(l-T))  CJX Oa)  HQ(Pa)  -JQ(Pa)  H^pa)  ] 

2it0 

B1  - - [-41)  J1(3a)  ^(aa)  +  (1+T1)  (paH^oa)  JQ(pa) 

+  aaJ^Pa)  HQ(aa))  -  aPa2JQ(Pa)  HQ(aa)  ] 


(33) 


n  *  2 

2|0ninpa 

A  - - - -  [J  (pa)H  (pa)  -  J  (pa)H  (pa)] 

n  An  n-1  n-1  n 

n 


B 

n 


[naaH  (aa)J(Pa)  +  npa  H  (aa)J  , (pa) 
n-1  n  n  n-i 


-  opa2  H  . (aa)J  r(pa)] 
n-i  n-i 


(34) 
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where 

6X  =  -4TlH1(aa)H1(pa)  +  (1+Tp  aaHgfaa)!^  (pa) 

+  (1+T|)  0aHo(0a)H1(aa)  -  a3a2H0(aa)HQ(3a) 

(35) 


A  =  naa  H  (aa)  H  (0a)  +  n0aH  .  (0a)H  (aa) 
n  n-l  n  n-1  n 

-  a0a2H  .  (aa)H  (0a) 
n-i  n-i 


(36) 


Substitution  of  Equations  (32)  through  (36)  into  Equations  (10)  and 
(11)  yields  the  solution  for  the  stresses  and  displacements  for  the 
rigid  inclusion  case. 

At  the  boundary,  where  r  =  a,  the  dimensionless  stresses 

ic  &  k 

T  >  and  t  .  are  simplified  to 

rr  00  r9  v 


k 

T 

rr 


(l-Tl)H1(aa) 


sin  6  +  ^ 
n=2 


in+*n  H  (aa) 
n 


sin  n0  ] 


(37) 


* 


T 


00 


k 

T 

rr 


(38) 


1 6 


2  2 

T*  .2  f  ,,lp  g.Z 

r0  rr  l 


2  2 


33a3H1(pa)  +  871  (  H0(pa) -pal^  (pa)) 


-  —  [  (1+T])H^  (aa)  -aaHQ(aa)  ]  cos  0 


-2 


n=2 


in+^(-nH  (aa)  4-  aaH  -  (aa)  ) 

n  n-l  A 

-  cos  n0r 


(39) 


It  is  to  be  noted  that  the  Poisson's  ratio  range  for  common  materials 

is  0  <  a  <  k;  therefore,  taa  is  always  less  than  t  .  The 

00  '  rr 

dimensionless  rigid  body  translation  and  rotation  of  the  inclusion 
are  given 

•  -ei^-[Hi(o*) -f  <40) 


*  16171 _ 

nPa[p2a2H1(pa)  +  8T)  (  ^  HQ(pa)  -H^pa))] 


(41) 


The  rotation  is  non-dimensional ized  with  respect  to  the  dimensional 
coefficient  of  the  rotation  associated  with  the  incident  wave,  i.e., 


n 


o 


2 


IT 


LONG  WAVELENGTH  LIMIT 

If  the  wavelength  of  the  incident  shear  wave  approaches  infinity 
(|3-*0) ,  the  solutions  obtained  for  the  cavity  and  the  rigid  inclusion 
should  reduce  to  the  Kirsch's  and  Goodier's  static  solutions.  Under 
pure  shear  Tq,  the  static  solution  for  the  cavity  case  at  the 
boundary  is 


—  =  T*  =  -4  sin  20 
tq  96 


(42) 


and  the  static  solution  for  the  rigid  inclusion  is 


rr  _  * 

tq  "  Trr 


2k 

k2+l 


sin  20 


T  2 

_r0  _  *  =  2k 

To  "  Tre  "  k2+l 


COS  20 


(43) 


*  ,i  2  .  * 

T00  =  (1  "  71  )  Trr 

k 


For  small  arguments,  the  Bessel  functions  are  approximated  by 
the  leading  terms  of  their  series.  Thus,  as  x  -»  0, 


Hn(x)  -4  1  +  i  (  -  )  (In  x-0. 11593) 
u  tt 


Hn(x)  -  (  f  )%  -i(n-l) |  )% 
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The  coefficients  of  sinusoidal  functions  for  in  the  cavity 

case  are  reduced  in  the  limit,  as  3a-*0  and  aa  =  ^  -♦  0,  to 

the  following  form  for 


n  =  1 


Q2  2 

~2~  HjCaa) 


(44) 


n  ^  2 


a  2  2 

n(n2-l-  ,a  )H  (aa) 

_ / _ n _ 


-  iTTti  (n2  -1-  4s"  )  (n-1) '  (  &  )n"2  /  4[(n-l)'.]2(n2+n-  ^  ) 


2  2  2  2  2 
+  is  -l)fl  a  £  (n_2)  .  f  +  2  (_n3+n_  Ls-  )  (n.2) !  (n-1) (1+  \  ) 

k 
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t 


- 


n  *  2 


V  0 


n  >  2 


(45) 


Therefore,  as  0a -*  0, 


Tee 


-4  sin  26 


which  is  identical  to  Equation  (42). 


(46) 


The  coefficients  for  the  sinusoidal  function  in  Equations  (37) 
and  (39)  for  the  rigid  inclusion  case  are  reduced  in  the  limit,  as 
aa  -»  0  and  |3a  -*  0,  to 


•Cl-TpHjfea) 


ind-TOBa 

871 


nH  (oa) 
n 

A 

n 


.  , 2  .  Ba  Nn-2 
ink  (  2  ) 

(n-2)  •  (1+k2) 


2  2 

B  a  (n-2)'. 
n(n-l) 


20 


ink 


2(l+k2) 


n  -  2 


n  >  2 


(47) 


2  2 
■IB  a 


2  2 

33a3H1(pa)  +  8T|  HQ(3a)  -  3a  H^Ba)] 


_  _n_  e22 

16T]  P 


(1+T|)  H^Caa)  -  aa  HQ(aa) 


-  in(l+Tl)Ba 

87] 


(48) 


-n  H  (aa)  +  aa  H  .  (aa) 
n _ n-i _ 


.  , 2  ,  Ba  >n-2 
-ink  (  *y  ) 

(n-2) 12  (1+k2) 


■ink 


2  (1+k2) 


n-2 


n  >  2 


Therefore , 


*  2  k 

T  -  *f—  sin  29 
09  kZ+l 
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*  2k 

T  —  cos  20 

r9  k  +1 

which  are  identical  to  Equation  (43). 

For  the  case  where  T)  =  0,  a  singularity  will  occur  in  the 
n  =  1  term  as  0*  0.  The  coefficients  in  this  case  are  of  the 

form 


H^aa)  k2 

L1  (l+k2)3a 

-H^  (aa)  -aa  (aa)  _^2 

A1  (l+k2)pa 


(49) 


Hence,  it  is  apparent  that  the  stresses  will  become  infinitely  large 
as  the  wave  number  approaches  zero  for  the  "fixed  in  space"  rigid 
inclusion. 


NUMERICAL  RESULTS  AND  DISCUSSIONS 

For  wavelengths  other  than  very  long  or  very  short,  the  stresses 
and  displacements  are  determined  by  the  series  given  in  Equations  (18) 
through  (20)  for  the  cavity  case  and  in  Equations  (37)  through  (41) 
for  the  rigid  inclusion  case.  Numerical  evaluations  were  carried  out 
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on  a  7090  IBM  computer  where  the  series  were  terminated  when  accuracy 

-4 

of  10  was  achieved. 

The  expressions  for  the  stresses  and  displacements  are  all  of  the 

form 


(R  +  iT) 


-i(Dt 


=  (R  + 


i2>* 


-i(tut-y) 


In  one  complete  cycle,  the  real  parts  R  yield  the  stresses  at  t  =  03 

T 

and  the  imaginary  parts  I  yield  the  stresses  at  t  =  where  T  is 

2  2  \ 

the  period  of  the  incident  wave.  The  absolute  values  (R  +  I  )  are 


the  maximum  stresses  which  occur  at  an  instant  between  t  =  0  and 


t  =  depending  on  the  phase  angle  y  =  tan 


■*(*)• 


Cavity 

In  this  case,  the  only  nonvanishing  stress  at  the  boundary  is 

* 

tqq,  therefore,  it  is  also  the  maximum  principle  stress.  Numerical  results 
of  stresses  and  displacements  vrere  obtained  for  the  values  of  a  «  0.15, 
0.25,  0.35  and  dimensionless  wave  numbers  in  the  range  0.1  £  pa  £  3.0. 

Fig.  2  illustrates  the  angular  distribution  of  HqqI  f°r  three 
wave  numbers,  pa  =  0.1,  1.0  and  1.5.  It  is  to  be  noted  that  at 

Pa  =  0.1,  the  stress  distribution  is  almost  identical  to  that  of  the 
static  case,  while  at  higher  wave  numbers  the  stress  distribution  is 
considerably  distorted  from  the  static  case.  Figs.  3  and  4  show  how 
stresses  at  two  points  on  the  cavity  (0  *  )  vary  with  frequency 
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2k 


and  the  Poisson's  ratio  of  the  medium.  In  this  case,  as  well  as 
results  shown  in  References  2,  3  and  4,  the  maximum  dynamic  stress  is  only 
10  to  20  per  cent  higher  than  the  static  values,  and  the  wave  number 
at  which  the  maximum  occurs  lies  between  0.20  and  0.50. 

The  angular  distributions  of  ju^J  and  (u^)  for  values  of 
pa  =  0.10,  1.0,  and  1.5  are  shown  in  Figs.  5  and  6.  At  pa  =  0.10  the 
displacements  are  again  the  same  as  the  static  case;  while  at  high 
frequency  there  appear  to  be  two  stationary  values  for  |u*|  in  the 


interval  of  0  £  9  £  tt,  and  three  stationary  values  for  .  |u^|  .  Both 
mraxima  for  |u*|  and  |u*|  tend  to  shift  toward  the  incident  side  as 
the  wave  number  increases . 
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Rigid  Inclusion 

Figs.  7  and  8  illustrate  the  angular  distributions  of  |t*  |  and 

"k 

| |  for  two  wave  numbers  (pa  =  0.1,  1.5)  and  7)  *  0.5.  At  pa  =  0.1 
(low  frequency) ,  the  distribution  is  nearly  symmetric  with  respect  to 
the  x  and  y  axes,  as  should  be  expected.  At  high  wave  numbers 
this  symmetry  does  not  exist,  and  the  maximum  tend"  to  shift  to  the  in¬ 
cident  side  of  the  inclusion. 


Fig.  7  Distribution  of  |Trr|for  Various  with  cr  =0.25 
and  7]  -  0.5  (Rigid  Inclusion) 


Figs.  9  through  12  show  the  variation  of  |r*r( ,  Jt*q|,  |©*|  ,  and 

|tl|  as  a  function  of  wave  number  for  various  values  of  7).  It  is 
7  P. 

apparent  that  the  density  ratio  7]  (— -  )  has  a  large  influence  on 
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Distribution  of  |  T  *^|  for  Various  fia.  with  a  =0.25  and  tj  =  0.5 
(Rigid  Inclusion) 
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the  stresses,  displacements,  and  rotation.  As  the  density  of  the  in¬ 
clusion  increases  the  maximum  value  of  the  stress,  displacement,  and 
rotation  also  increases.  This  increase  becomes  appreciably  greater 
for  density  ratios  below  0.5.  In  the  case  of  T)  =  0.25,  for  example, 
the  maximum  dynamic  stress  for  lTrrl  an<*  |  Tr q I  * n  Figs  •  9  and  10  are, 
respectively,  37  and  105  per  cent  greater  than  those  under  the  static 
condition . 


II  „  .  7T  37T 

Fig.  9  |  Trrl  vs  crt  "7  f°r  Various  rj 

with  cr  =  0.25  (Rigid  Inclusion) 


It  is  also  to  be  noted  that  there  is  an  apparent  relationship 
indicated  between  the  maximum  values  of  stresses  and  those  of  both  the 
rigid  body  motions.  Previous  studies^’  ^  have  shown  that  there 
exists  only  one  maximum  for  stresses  as  a  function  of  frequency,  which 
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is  also  the  case  for  lTrrl  Fig.  9.  However,  Fig.  10  shows  for 

I  *  I 

|t  q[  at  9  =  tt  that  there  are  two  distinct,  stationary  values  for 
Tj  =  0.10  at  two  different  wave  numbers.  By  comparing  these  wave 

i  *  i 

numbers  where  the  peaks  occur  with  those  where  the  maximum  of  [U  J 

and  J©  |  occur,  it  is  apparent  that  the  first  maximum  of  [t^) 

i*i 

occurs  very  close  to  the  wave  number  where  |U^|  is  also  a  maximum. 
Furthermore,  the  wave  number  corresponding  to  the  second  maximum  of 
|  |  is  very  close  to  the  wave  number  at  which  j©  |  is  a  maximum. 

These  observations  certainly  indicate  the  presence  of  coupling 


between  the  rigid  body  motions  and  the  stress  at  the  boundary. 


Fig.  11  ©*  vs  /3a  for  Various  17  (Rigid  Inclusion) 
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Fig.  12  |  U*|vs  P a  for  Various  7]  with  cr  =0.25 

(Rigid  Inclusion) 


From  Figs  .  11  and  12  it  is  shown  that  the  denser  inclusions 
receive  the  greater  rigid  body  motion  for  pa  approximately  less  than 

i  ★ .  i^i 

1.5  for  I®  I  and  1.0  for  |U^|  ,  but  the  converse  is  true  as  the  wave 
number  increases  beyond  these  values  . 
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